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ABSTRACT 

As part of a long term monitoring campaign of Mrk 335, deep XMM-Newton observations 
catch the narrow-Hne Seyfert 1 galaxy (NLSl) in a complex, intermediate flux interval as the 
active galaxy is transiting from low- to high-flux. Other works on these same data examined 
the general behaviour of the NLS 1 (Grupe et al. ) and the conditions of its warm absorber 
(Longinotti et al. ). The analysis presented here demonstrates the X-ray continuum and timing 
properties can be described in a self-consistent manner adopting a blurred reflection model 
with no need to invoke partial covering. The rapid spectral variability appears to be driven by 
changes in the shape of the primary emitter that is illuminating the inner accretion disc around 
a rapidly spinning black hole (a > 0.7). While light bending is certainly prominent, the rather 
constant emissivity profile and break radius obtained in our spectral fitting suggest that the 
blurring parameters are not changing as would be expected if the primary source is varying its 
distance from the disc. Instead changes could be intrinsic to the power law component. One 
possibility is that material in an unresolved jet above the disc falls to combine with material 
at the base of the jet producing the changes in the primary emitter (spectral slope and flux) 
without changing its distance from the disc. 

Key words: galaxies: active - galaxies: nuclei - galaxies: individual: Mrk 335 - X-ray: 
galaxies 



1 INTRODUCTION 

The narrow-line Seyfert 1 galaxy (NLSl), Mrk 335 {z — 
0.025785), has always been the subject of intense scrutiny. As one 
of the brighter active galactic nuclei (AGNs) in the X-ray sky it has 
been observed by most major missions so its long term behaviour is 
relatively well established (Grupe et al. 2008). Its X-ray continuum 
has normally been described with a power law and soft-excess (e.g. 
Pounds et al. 1987). More sensitive observations revealed a narrow 
emission line around 6.4 keV along with a broader emission feature 
in the Fe Ka region. The combination of the soft excess and broad 
emission feature is indicative of reflection from the inner region of 
an ionised disc (e.g. Ross & Fabian 2005) and such models have 
been successful in fitting the spectra of Mrk 335 (e.g. Ballantyne, 
Iwasawa & Fabian 2001; Crummy et al. 2006; O'Neill et al. 2007; 
Longinotti et al. 2007; Larsson et al. 2008). 



A 2007 Swift snap-shot of Mrk 335 found the NLSl in an 
unprecedented low state where the X-ray flux had diminished to 
about 1/30 of the previously lowest observed flux state (Grupe et 
al. 2007). The Swift observation showed considerable curvature in 
the X-ray spectrum, and simultaneous optical data found that the 
AGN was in an X-ray weak state. Both factors are indicative of ei- 
ther absorption or reflection dominating the spectrum (Gallo 2006). 
The Swift observation prompted a 20 ks XMM-Newton Target of 
Opportunity observation that revealed low-energy emission lines 
from a photoionised gas likely originating far from the inner ac- 
cretion disc and close to the broad line region (Grupe et al. 2008; 
Longinotti et al. 2008). During this low state the continuum could 
be described equally well as arising from either partial covering or 
blurred reflection (Grupe et al. 2008). 

Since 2007 Mrk 335 has undergone continuous monitoring 
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with Swift. The AGN has exhibited significant UV and X-ray vari- 
ability while remaining in a moderately low X-ray flux state. A 
deep 200 ks pointed XMM-Newton observation in 2009 attempted 
to catch the AGN in a minimtmi in order to study the properties 
of the photoionised gas observed in the 2007 low state; however, 
the AGN was instead caught in a transition as it brightened to an 
intermediate flux level (Grupe et al. 2012, hereafter G12). Inter- 
estingly, while the emitting gas is less evident in this intermediate 
state due to the higher continuum flux compared to the 2007 low 
state, a warm absorber is now clearly detected (G12; Longinotti et 
al. 2012, hereafter L12). Based on their analysis of the RGS data, 
LI 2 describe the warm absorber as a combination of three different 
zones with distinct temperatures and ionisation parameters. How- 
ever, the RGS analysis shows that the ionized gas does not produce 
the variability observed during this intermediate flux state. The cur- 
vature in the X-ray spectrum again could be described by partial 
covering prompting the suggestion that clouds of absorbing mate- 
rial have been entering the line of sight since 2007 and variations in 
covering fraction are responsible for the recent low state and vari- 
ability exhibited by Mrk 335 (G12). Grupe et al. (2008) point to 
Tanaka et al. (2005) for a partial covering explanation of Mrk 335 
in the high-flux state when the covering fraction would have been 
low. 

If the recent diminished flux is due to increased photoelectric 
absorption from additional partial covering, then it should be ac- 
companied by increased emission in the narrow Fe Ka line (e.g. 
Reynolds et al. 2009). G12 find that within uncertainties the Fe Ka 
flux in Mrk 335 is identical in the intermediate- and low-flux states, 
both of which are comparable (perhaps even lower) than the Fe Ka 
flux reported by Grupe et al. (2008) in the high-flux state. In addi- 
tion, one may expect to find markers of absorption in the UV. There 
is some indication of absorption in a recent HST-COS spectrum of 
Mrk 335 (L12). However if and how this is linked to any absorp- 
tion in the X-ray is difficult to establish since the UV and X-ray 
observations were not simultaneous. 

In this work we test the blurred reflection interpretation for 
Mrk 335 when the AGN is in an intermediate low-flux state and is 
significantly modified by a warm absorber, to determine if partial 
covering is necessary to describe the behaviour of the source in 
recent years. The observations and data reduction are described in 
the following section. Spectral modelling is explored in Section 3 
and the temporal analysis is described in Section 4. We discuss and 
summarize our results in Section 5 and 6, respectively. 



2 OBSERVATIONS AND DATA REDUCTION 

Mrk 335 has been observed with XMM-Newton (Jansen et al. 2001) 
on five occasions between 2000 and 2009. In this work we report 
on the long 200 ks observation spread over two consecutive or- 
bits (revolutions 1741 and 1742) in 2009. During the observations 
the EPIC pn detector (Striider et al. 2001) was operated in large- 
window mode and with the thin filter in place. The MOS (MOS 1 
and MOS2; Turner et al. 2001) cameras were operated in full-frame 
mode and with the medium filter. The Reflection Grating Spectrom- 
eters (RGSl and RGS2; den Herder et al. 2001) also collected data 
during these observations, as did the Optical Monitor (CM; Mason 
et al. 2001). The OM data have been presented in G12 and will not 
be reanalyzed in this work. These data will be discussed as needed. 
The RGS analysis is the subject of an in-depth study of the warm 
absorber (L12). We will incorporate the results of L12 in this work, 
but we will not reanalysis the RGS data here. 



The XMM-Newton Observation Data Files (ODFs) from each 
observations were processed to produce calibrated event lists using 
the XMM-Newton Science Analysis System (SAS vl 1 . . 0). Un- 
wanted hot, dead, or flickering pixels were removed as were events 
due to electronic noise. Event energies were corrected for charge- 
transfer inefficiencies. Light curves were extracted from these event 
lists to search for periods of high backgroimd flaring and such peri- 
ods have been neglected. Pile-up was negligible during the observa- 
tions. The background photons were extracted from an off- source 
region on the same CCD. Single and double events were selected 
for the pn detector, and single-quadruple events were selected for 
the MOS. EPIC response matrices were generated using the SAS 
tasks ARFGEN and RMFGEN. The MOS and pn data at each epoch 
were compared for consistency and determined to be in agreement 
within known uncertainties (Guainazzi et al. 2010) as was the case 
in G12. For simplicity we will only present the pn data at each 
epoch. The total amount of good exposure for the pn during revo- 
lution 1741 and 1742 is 96.4, and 70.0 ks, respectively. 

The XMM-Newton X-ray spectra were grouped such that each 
bin contained at least 30 counts. Spectral fitting was performed us- 
ing XSPEC vl 2 . 7 . 1 (Amaud 1996). AU parameters are reported 
in the rest frame of the source unless specified otherwise. Figures 
remain in the observed frame and are binned for display purposes 
only. Unless stated otherwise, black, open symbols identify data 
from revolution 1741; and filled, red squares represent data from 
orbit 1742. The quoted errors on the model parameters correspond 
to a 90% confidence level for one interesting parameter (i.e. a Ax^ 
= 2.7 criterion). A value for the Galactic column density toward 
Mrk 335 of 3.56 x 10^° cm^^ (Kalberla et al. 2005) is adopted in 
all of the spectral fits and abimdances are from Anders & Grevesse 
(1989). Limiinosities are calculated using a Hubble constant of 
Ho=70 km Mpc~^ and a standard flat cosmology with S7m 
= 0.3 and = 0.7. 



3 SPECTRAL MODELING 

3.1 The high-energy spectral region 

The 2.5—10 keV band in Mrk 335 typically displays curvature both 
in the high- and low-flux states that can be attributed to blurred re- 
flection (e.g. O'NeiU et al. 2007; Larsson et al. 2008; Grupe et al. 
2008) or partial covering (e.g. Grupe et al. 2008; G12). Narrow 
features have also been detected. A narrow ~ 6.4 keV emission 
line arising from Fe Ka is regularly observed, while narrow emis- 
sion features from a highly ionized plasma (Fe XXVI Ly-a, O'Neill 
et al. 2007) or absorption features attributed to relativistic inflows 
(Longinotti et al. 2007) are transient. 

The 2.5 — lOkcV spectrum from each epoch was fitted with 
a power law and a second component to mimic distant reflection. 
For the distant reflector we consider two models, pexmon (Nandra 
et al. 2007) and ref lionx (Ross & Fabian 2005). Both models 
self-consistently include Fe Ka, the Compton shoulder, and Comp- 
ton hump. In addition, pexmon includes Fe K/3 and Ni Ka emis- 
sion. However, pexmon is for reflection off neutral material and 
does not consider the features that arise below the iron-K band. For 
ref lionx the material is also considered nearly neutral as the 
ionisation parameter (^ = ArcFx/n, where is the incident flux 
and n is the density of the gas) was set to Icrgcms^^ (the low- 
est allowed by the model). Iron abundance is also fixed to the solar 
value. While there is evidence that the inner accretion disc can show 
super-solar iron abundances (e.g. Fabian et al. 2009; Reynolds et al. 
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Figure 1. The 2.5 — lOkeV spectrum is fitted at each epoch with shghtly 
different models for the distant reflector. Top panel: Residuals (data/model) 
remaining in the 3.5 — 7.5 keV after fitting the spectra with a power law 
plus pexmon model. Lower panel: Same as above, but using ref lionx to 
model the narrow Fe Ka emission Hne. In both figures the black open circles 
con'espond to the data from revolution 1741 and the filled, red squares are 
data from revolution 1742. This colour/symbol convention is used through- 
out the paper unless stated otherwise. 



2012) it is not obvious that this characteristic extends to the distant 
reflector. However, we note that Walton et al. (2012a) was able to 
model the Suzaku spectrum of Mrk 335 when they allowed the dis- 
tant reflector (torus) and the inner accretion disc to share common 
abundances. 

Statistically the models fit reasonably well (x^/dof = 
1.03/1507 and 1.01/1508 for the pexmon and reflionx, re- 
spectively), but both models show excess in the residuals between 
3.5 — 7.5 keV (Fig.[TJ. G12 interpreted this excess as evidence for 
partial covering, we examine this in terms of blurred, ionised re- 
flection from the inner accretion disc. Adding a single Laor pro- 
file to both models improved the fits (xt < 1-0) and residu- 
als significantly in both cases. The equivalent width of the line is 
EW ~ 400 eV in both models. The fitted line parameters are the 
energy (E ~ 6.4 keV), disc emissivity index (g ~ 4 — 8), inner 
disc radius (Ri„ ~ 1.3 — 5rg, [rg = GM/c^]), and disc inclina- 
tion (z ~ 30 — 65°). We do not examine these parameters in detail 
at this point as they will serve as initial values for more realistic 
models in the remaining analysis. 



3.2 The mean 0.3 12 keV spectrum 

As seen in Fig. [2] the average spectrum undergoes significant 
changes between revolution 1741 and 1742. The difference be- 
tween the two spectra (i.e. rev. 1742 - rev. 1741) fitted with an 
absorbed power law (F « 2.5) between 2 — lOkeV and extrapo- 
lated to lower energies is shown in the lower panel of Fig. |2l The 
changes are substantial below ~ 1.5 keV, but there is also indica- 
tion for variability above 2 keV. For this reason we will treat the 
data from each epoch separately when examining the spectra. 

There is a significant warm absorber present at both epochs, 
which is clearly seen in the RGS spectra, and is the subject of other 
work (LI 2). This is the first time that a warm absorber has been un- 
ambiguously observed in Mrk 335. L12 identified three distinct ion- 




Observed energy (keV) 

Figure 2. Top panel: The average specfium in revolution 1741 (black cir- 
cles) and revolution 1742 (red squares). Lower panel: The ratio of the dif- 
ference spectrum (rev. 1742 - rev. 1741) fitted with an absorbed power law 
(F fa 2.5) between 2 — 10 keV and extrapolated to lower energies. As seen 
in the top panel the vaiiability is dominate below 1.5 keV. 



ized components, all outflowing at high velocity (~ 5000 km ). 
However, since L12 determined that the warm absorber itself is not 
variable during these observations, it cannot be the warm absorber 
that is responsible for the changes seen in Fig. |2] To remain con- 
sistent with the RGS analysis, we constructed XSTAR (Kallman & 
Bautista 2001) grids that are comparable with the three warm ab- 
sorbers found in L12. These components are included in our spec- 
tral fitting and are linked between the epochs. 

The neutral reflector originating from distant matter is mod- 
elled with reflionx rather than pexmon since it offers a more 
complete treatment of the reflection at low energies. This compo- 
nent should not be confused with the reflection from the inner disc 
as it is not subject to any blurring effects. The ionisation parameter 
of this distant reflector is fixed at ^ — 1.0 erg cm ^ (i.e. the low- 
est value permitted by the model) and the iron abundance is fixed 
to solar. The photon index of the ionising power law is fixed to the 
canonical value of F = 1.9 (Nandra & Pounds 1994). It is not clear 
if this assumption is accurate as the photon index is clearly variable 
in Mrk 335 (e.g. 1.5 < F < 2.3 have been reported). It is also 
not obvious that the torus would be illuminated by the same power 
law that is seen by the observer since the line-of-sight through the 
corona, and hence the optical depth is likely different from each 
perspective. The normalisation of this component is linked at both 
epochs. 

In the 2007 low flux state observation, Grupe et al. (2008) 
and Longinotti et al. (2008) identified narrow emission features at- 
tributed to a distant photoionised plasma. The CCD spectra were 
fitted with a MEKAL model to account for this emission (Grupe et al. 
2008), but the fitting was rather ad hoc and did not account for some 
of the strongest features in the spectrum (e.g. a 0.89 keV emission 
feature and the narrow Fe Ka emission line). Such emission is not 
prominent in the 2009 RGS data (LI 2) and this could simply be a 
matter of the higher continuum flux overwhelming the hot, diffuse 
plasma. We do not specifically consider the reasonable assumption 
that a constant MEKAL component (e.g. fixed to the 2007 values) 
is present in these 2009 data, but we note that the distant reflector 
modelled by the unblurred reflionx component can similarly 
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Figure 3. Top panel: The average spectrum in revolution 1741 (black cir- 
cles) and revolution 1742 (red squai'es) over the 0.3 — 12keV band fitted 
with the best-fit model (Table[T). The dotted lines are the blurred reflector 
and the dashed curves are the intrinsic power law. The distant reflector is 
common to both spectra and is shown in blue (dashed-dotted curve). Lower 
panel: The data-to-model ratio resulting from the best-fit model. 



describe this emission. The 0.5 — 10 keV flux of the MEKAL com- 
ponent in the low-flux spectrum was ~ 2.2 x 10~^^ erg cm^'^ , 
comparable to the flux of the distant reflector used in these models 
(see TablelTJ. 

After fitting the three-phase warm absorber to the RGS data, 
L12 attributed the remaining residuals between about 0.49 — 
0.65 keV to intrinsically narrow emission corresponding to O VIII 
and N VII Ly-a, as well as the O VII triplet. Adding Gaussian pro- 
files with fixed energies and widths, but free normalisations, to the 
CCD data did not result in any detection of these features. For each 
line, the measured equivalent width was EW < 1 eV. Including 
the Gaussian features with all parameters fixed to the RGS mea- 
sured values generated an entirely comparable fit to the one pre- 
sented in Table [T] 

The broadband continuum is described with a power law plus 
reflection from the accretion disc. It is then modified by a three- 
zone warm absorber and a distant neutral absorber. The reflection 
component is blurred for dynamical effects in the disc close to the 
black hole using kdblur2 with a broken power law emissivity 
profile. The objective was to obtain a self-consistent model that de- 
scribed both observations. Despite the significant variability the ob- 
servations were separated by only ~ 42 ks, consequently we were 
selective in what we allowed to vary between the observations. In 
addition to the distant reflector, the inclination (i), iron abundance 
(Ape), inner (Rin) and outer {Rout) disc radius were linked be- 
tween the observations. Initial fits were conducted with the inner 
emissivity profile (qi„) and the break radius (Rt) free to vary at 
each epoch. The outer emissivity profile (qout) was fixed to 3 as 
would be expected in the case for illumination of the disc from a 
distant source. Finding no significant differences in the values of 
Rb and qi„ during revolution 1741 and 1742, the parameters were 
also linked in subsequent fits. The best-fitting model reported here 
allow only the photon index and normalisation of the power law 
component; and the ionisation and normalisation of the reflector to 
vary at each epoch. 

The model adequately fits the data (x^/dof = 1.06/2379). 



Figure 4. The best-fit model described in this section is extrapolated to high 
energies and compared (i.e. not fitted) to the Swift BAT data . The residuals 
are shown in terms of sigma. 




Spin parameter (a*) 

Figure 5. The quality of the fit is tested against the black hole spin pa- 
rameter showing that extremely low spins are discriminated against. The 
minimum m Ax^ is at a 0.9. The green, dashed line marks the 90 per 
cent confidence level. 

The model parameters and fit are shown in Table [T] and Fig. [3] re- 
spectively. The reflection model makes specific predictions of the 
spectral appearance above lOkeV. While simultaneous observa- 
tions above lOkcV do not exist, the average Swift BAT spectrum 
is available for comparison (Fig.|4l(. The Swift data agree well with 
blurred reflection spectral model. 

The differences between the two epoch require changes in the 
normalisations as well as the shape of the reflection and power law 
components. In contrast to typical observations the reflection com- 
ponent appears more dominant during revolution 1742 when the 
AGN is in a higher flux state. 

3.3 Black hole spin 

In accordance with the small inner radius found in the fit (Table [TJ 
a high black hole spin value is expected in Mrk 335. In order to es- 
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Table 1. The best-fit model parameters for the multi-epoch spectral modelling of Mrk 335. The model components and model parameters are listed in Columns 
1, and 2, respectively. Columns 3 and 4 list the parameter values during revolution 1741 and 1742, respectively. The inner (Rin) and outer (Rout ) disc radius 
are given in units of gravitational radii (1 rg = GMjc?). The reflection fraction (7?.) is approximated as the ratio of the reflected flux over the power law 
flux in the 20 ~ 50 kcV band. Values that are linked between epochs appear only in column 3. The superscript / identifies parameters that are fixed and the 
superscript p indicates values that are pegged at a limit. Fluxes are corrected for Galactic absorption and are reported in units of erg cm^"^ s^^ . 
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timate the black hole spin we conservatively concentrate efforts on 
the high-energy spectrum above 2.5 keV. When fitting the broad- 
band spectrum the spin value is often highly dependent on the inter- 
pretation and modeling of the smooth soft excess where the statis- 
tics are highest, and as we know from LI 2, this region is signifi- 
cantly modified by the three-zone warm absorber. 

We begin with the power law plus distant reflector model 
described in Section 13.11 To this we add a reflector that is 
blurred using kerrconv (Brenneman & Reynolds 2006; instead 
of kdblur2), which includes black hole spin as a model param- 
eter. As in the broadband spectral analysis, all model parameters 
are tied except for the power law photon index and normalisation, 
and the blurred reflector ionisation and normalisation. A good fit 
to the 2.5 - lOkeV data is found (xl/doi = 0.96/1307) with 
all model parameters in agreement with the broadband results in 
Table[T](albeit with larger uncertainties). 

In order to estimate the uncertainty on the measured spin we 
examine the variation in the quality of the fit associated with the 
changes in the spin parameter (Fig. O. Relatively high spin val- 
ues are preferred by the fit (0.76 < a < 0.998) at 90 per cent 
confidence. Very low spin values (a < 0.2) are rejected at the So- 
level. The spin parameter measured from these data of Mrk 335 
in a warm absorber dominated, intermediate flux state is in agree- 
ment with that measured by Walton et al. (2012a) using Suzakii 
XIS and PIN spectra of Mrk 335 in an unabsorbed, high-flux state 



(a = 0.8±0.1), and commensurate with the small inner disc radius 
found by Crummy et al. (2006). 



4 X-RAY VARIABILITY 

G12 reported significant flux and spectral variability over the inter- 
mediate flux level observations examined here (see their Fig. 5). We 
attempt to describe the variability in Mrk 335 based on the blurred 
reflection model presented in Section [T2l As noted in Section [J!2l 
variability is seen in power law normalisation and shape as well as 
in the ionisation and flux of the reflector. Here we examine how the 
model can account for the spectral variability as a function of time 
and flux level. 



4.1 Flux resolved spectroscopy 

Spectra are created in four flux bins during revolution 1741 and 
three flux bins during revolution 1742, as defined in Fig. [6] The 
high-, medium-, and low-flux states are defined identically in 1741 
and 1742. The extra-low-flux state (i.e. xlow) is unique to 1741 as 
the AGN brightness does not dip to this level during 1742. In Fig.|7] 
we plot the spectra from each flux state against a power law (F = 
2) absorbed by Galactic column density. Despite the comparable 
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Figure 6. The 0.2 — 12 keV pn light curve in 500 s bins from revolution 
1741 (left) and 1742 (right). Zero second marks the start of each observa- 
tion. The dashed lines delineate the various flux levels examined in Sec- 
tion |4]T] The high, medium, and low state are defined with count rates of 
> 6, 4.5 — 6.0, and 3 — 4.5 count . Revolution 1741 also is separated 
into an extra low state (< 3 count s~ ^ ). 
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Figure 8. The ratio (data / model) from fitting the various flux states defined 
in Fig.|6]with a model in which the reflection and power law component are 
allowed to vary (see text for details). Data from 1741 and 1742 are shown 
in the top and bottom panel, respectively. Black, red, and green data points 
mark the low-, medium-, and high-flux states. Blue data points in the top 
panel identify the exfia-low-flux state during 1741. 
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Figure 7. Spectra extracted from the flux levels defined in Fig.|6]are com- 
pared to an absorbed power law with F = 2. Despite the comparable flux 
there exist differences between the flux-resolved spectra from the two ob- 
servations. 



brightness the flux-resolved spectra differ notably between 1741 
and 1742. 

The various flux states during each orbit were fitted simul- 
taneously using the best-fit model from Section [X2l as the starting 
condition. The warm absorber and distant reflector parameters were 
fixed to their average values. Initially, the blurring parameters qin 
and Rh were free, but since neither demonstrated noticeable vari- 
ability they were fixed to their average values as well. The only 
free parameters were the photon index of the power law contin- 
uum, ionisation parameter of the reflector, and the normalisation of 
both components. 

The spectral variations are well fitted by allowing only these 
parameters to vary (xi /dof = 1.02/3525 for revolution 1741 and 



xJ/dof = 1.01/2502 for revolution 1742). The residuals from 
the fits are shown in Fig.[8] Some trends resulting from fitting the 
flux-resolved spectra are shown in Fig. |9] As found when mod- 
elling the average spectra, the ionisation parameter of the blurred 
reflector is distinctly higher during revolution 1742. This is seen in 
the flux-resolved spectra as well; however, there does not appear 
to be a gradual change with increasing flux that might be expected 
(Fig|9^). We note the telescope does not collect data for ^ 40 ks 
between revolution 1741 and 1742, which makes it difficult to draw 
strong conclusions as to when the change in the ionisation parame- 
ter occurs. 

Other parameters do appear to exhibit predicable behaviour. 
The power law slope is seen to steepen with increasing flux 
(Figjgjj). The power law flux and reflected flux are also correlated, 
which would occur if the primary emitter is brightening intrinsi- 
cally and illuminating the disc more (Fig[9]:). The reflection frac- 
tion {TZ), which is taken to be the ratio between the reflected and 
power law flux between 20 — 50 keV, is also correlated with the 
reflected flux (Fig|9jl). We note the relationship between the two 
parameters is different at each epoch, which could be due to the 
simultaneous changes in the power law photon indices. 

4.2 Time resolved spectroscopy 

In order to investigate the model variability over time, spectra were 
created using 10 ks intervals of data resulting in 12 spectra in 1741 
and 7 spectra in 1742. Each spectrum was modelled by allowing 
only the power law normalisation and slope to vary as well as the 
reflector normalisation and ionisation. Fig. 1 101 and 1 1 1 1 portray our 
findings. 

The average ionisation of the reflector is different at each 
epochs, but as we see in the top panels of Fig. [TO] the ionisation 
parameter does not vary with the short timescale changes in the re- 
flected flux (lower panels of FigllOll. Similar behavior was seen in 
the flux-resolved spectra (Fig|9^). 

The flux of the reflector (lower panels of FigllOb is variable 
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Figure 9. The variability of various parameters in different flux states during 1741 (black open circles) and 1742 (red filled squares), (a) The ionisation 
parameter of the reflector as a function of the reflected flux (F^ ) in the 0.3 ~ 10 keV band, (b) The power law photon index as a function of power law flux 
(Fy) in the 0.3 — 10 keV band, (c) The reflected flux is plotted against power law flux, (d) The reflection fraction, between 20 — 50 keV, is compared to the 
reflected flux in the same band. 



over the course of the observations, as is the flux and the photon in- 
dex of the power law component (Figll It. Of interest, is that during 
the first 70 ks of 1741 the flux of the reflection component is ap- 
proximately constant (Fig.[TOl lower left panel) as is the shape and 
flux of the power law component (Fig [TT] left). In general, the re- 
flected flux appears to be well correlated with changes in the power 
law component. The specific correlation between the photon index 
and the reflected flux is shown in Fig. [12] and exists in both the 
time-resolved and flux-resolved spectra. Since the reflection frac- 
tion (TZ) increases with reflected flux (Fig|9jl), the observed trend 
in Fig[T2]between the shape of the intrinsic power law illuminating 
the accretion disc and the reflected flux is expected (see e.g. Figure 
4 ofRoss& Fabian 2005). 

The variability behaviour also describes the appearance of the 
fractional variability (Fvar) spectrum of Mrk 335 at each epoch 
(Fig. |13t . The Fvar in each energy band is calculated with light 
curves in 500 s bins following Edelson et al. (2002) and uncertain- 
ties are estimated following Ponti et al. (2004). There is clearly 
a difference below ^ 1 keV in the Fvar spectra between the two 
epochs with 1741 showing larger variations at lower energies. The 



Fvar spectrum is also calculated for only the first 70 ks of revolu- 
tion 1741 (shaded, black circles), showing much smaller variations 
over this reduced period. Each F^ai spectrum can be reasonably de- 
picted by the variations found in Fig. [lO] and [TT] The first 70 ks of 

1741 are described by variations in only the power law normalisa- 
tion, while to describe the average Fvar in 1741 (black open circles) 
variability in the remaining ~ 50 ks requires additional changes in 
the power law slope and reflector normalisation. The variability in 

1742 also requires changes in the power law shape and normalisa- 
tion, and in the reflector normalisation. 



5 DISCUSSION 

Blurred reflection is expected when the inner accretion disc extends 
down toward the radius of marginal stability around a black hole. 
Many of the concerns raised over blurred reflection are now com- 
ing into focus with recent X-ray data. The discovery of broadened 
Fe La and short lags (e.g. Fabian et al. 2009; Zoghbi et al. 2010) 
suggest the soft excess in AGN is due to reflection and that it arises 
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Figure 10. The variability of tlie reflection component is tracked over tlie entire observation by fitting 10 ks resolved spectra during revolution 1741 (left 
panel) and 1742 (right panel). The flux (lower panels) and ionisation parameter (upper panels) are compared to the average values (solid lines). The ionisation 
parameter appears constant during both observations, but the flux is variable. The reflected flux during revolution 1741 (lower left panel) is compared to the 
average during the first 70 ks of the observation (dashed line). The flux appears constant during this segment of the observation. 
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Figure 11. The variability of the power law component is tracked over the entire observation by fitting 10 ks resolved spectra during revolution 1741 (left 
panel) and 1742 (right panel). The flux (lower panels) and photon index (upper panels) are compared to the average values (solid lines) and both parameters 
appear variable. The photon index during revolution 1741 (upper left panel) is compared to the average during the first 70 ks of the observation (dashed fine). 
The photon index appears constant during this segment of the observation. 



from reprocessing of the primary continuum. The broad Fe Ka line 
is also seen in black hole binaries (BHB) (e.g. Miller et al. 2009) 
showing that reflection is ubiquitous in black hole systems. The 
similarities in the Fe Ka profile seen in AGN and BHB (Walton 
et al. 2012b) argues strongly in favour of the reflection interpreta- 
tion. Moreover, the evidence for blurred reflection is not limited to 
unique systems or specific lines-of-sight. Evidence of reflection is 
seen in most systems where there is a direct view of the central re- 
gions (e.g. Crummy et al. 2005; Zoghbi et al. 2012; Walton et al. 
2012a; De Marco et al. 201 1). 

The XMM-Newton observations of Mrk 335 analyzed here 
provide a view of the AGN that is distinct from previous obser- 
vations. Firstly, by the presence of a complex warm absorber (ana- 



lyzed by LI 2), and secondly through the significant and rapid spec- 
tral variability. Over the course of the 200 ks observation the source 
becomes notably softer as it brightens. When examining the light 
curve from the long-term Swift monitoring campaign (G12), we 
find that these XMM-Newton observations catch Mrk 335 in rapid 
transition from a low-flux to a high-flux state (G12). The X-ray 
continuum is well described with a primary power law emitter and 
blurred reflection from inner region of the accretion disc around 
a rapidly spinning black hole. The blurred reflection model also 
accurately describes the variability assuming that the continuum 
components are allowed to vary. Specifically, the photon index and 
power law normalisation as well as the reflector ionisation and nor- 
malisation are changing. 
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Figure 13. The fractional variability as of function of energy is calculated 
for observation 1741 (black open circles) and 1742 (red filled squares), as 
well as for the first 70 ks of observations 1741 (black, shaded circles). The 
models described in Section |4!2l are overplotted and show reasonable agree- 
ment with the data. 



spect to the disc (e.g. Reynolds & Fabian 1997; Beloborodov 1999) 
could be invoked to drive the changes in F. This interpretation has 
been advocated for other AGN and NLSls (e.g. Gallo et al. 2007, 
2011; Miniutti et al. 2010). However, while light bending is cer- 
tainly prominent, the rather constant emissivity profile and break 
radius obtained in our spectral fitting of Mrk 335 suggest that the 
blurring parameters are not changing as would be expected if the 
primary source is varying its distance from the disc. 

Fabian et al. (2012) suggest a double power law continuum 
to describe the low state X-ray spectrum of Cyg X-1. There they 
suggest the steeper component originates at the base of the jet 
and is responsible for the reflection with the harder component lo- 
cated farther up the jet. This scenario could perhaps be expanded 
to describe the observations in Mrk 335. One could imagine mul- 
tiple jet components coming together at the base of a jet located 
close to the disc, for example a jet "blob" that does not obtain the 
escape velocity and falls back downstream toward the base (e.g. 
Ghisellini et al. 2004). The coming together of two components 
of differing densities could produce the changes seen in the power 
law and ultimately in the reflection. Interestingly, a double power 
law continuum was suggested for the blurred reflection interpre- 
tation of the very low flux state observation of Mrk 335 (Grupe 
et al. 2008). Mrk 335 is certainly radio-quiet, but not radio silent 
(R = LviJjGH z) / Li,{B) ~ 0.25) so the possibility of some jet 
emission does exist. 



6 CONCLUSIONS 

The most recent XMM-Newton observations of Mrk 335 catch the 
AGN in a complex, intermediate flux state. Both the X-ray con- 
tinuum and timing properties of the NLSl can be described in a 
self-consistent manner adopting a blurred reflection model with no 
need to invoke partial covering. The rapid spectral variability ap- 
pear to be driven by changes in the spectral shape of the primary 
emitter that is illuminating the inner accretion disc. 
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Though the behaviour is complex, the variability in the pa- 
rameters is largely consistent with what is expected. The overall 
changes appear to be driven by the pivoting of the power law com- 
ponent (i.e. the primary emitter). The power law flux increases as 
the spectrum steepen, in line with expectations (e.g. Haardt et al. 
1997). Assuming the distance between the primary emitter and the 
disc does not change, the flux and ionisation of the disc will in- 
crease accordingly, as is seen. 

The changes in power law photon index occur very rapidly. 
During the first 70 ks of the observation Mrk 335 did not exhibit 
any notable spectral variability. The apparent change in the photon 
index, specifically steepening from F « 1.83 to F « 1.98, occurs 
over ~ 50 ks in the rest frame. A scenario in which the primary 
emitter could be moving at mildly relativistic velocities with re- 
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